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It is well known that the Casimir force between two half-spaces is dictated by their electromagnetic 
properties. In particular, when one of the half-spaces is mainly metallic or dielectric and the other is 
mainly magnetic, it is possible to show that the force is repulsive. This has attracted lots of interest 
towards the study of metamterials (MMs) in the context of Casimir effect, as their magnetic activity 
might help bring the idea of Casimir repulsion from the theoretical realm to experimental verification. 
In this paper we investigate the possibility of repulsion when the MM magnetic permeability is 
given not by a Drude-Lorentz behavior, but by a model put forward by Pendry et al. [20I ]. After 
introducing the model and deriving the necessary formulas, we show that it is impossible to achieve 
repulsion with such a model and present a qualitative discussion of why this is so. 
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. The last decade has witnessed an increased interest in Casimir physics [l|, Q thanks to improved precision measure- 
ments Q of the force between material objects separated by micron and sub-micron gaps. In a few words, the Casimir 
force may be thought of as a consequence of changing the vacuum fluctuations spectrum through the insertion of 
material boundaries. While the Casimir force offers new possibilities for nanotechnology, such as actuation mediated 
by the quantum vacuum, it also presents some challenges, as micro and nanoelectromechanical systems (MEMS and 
> ■ NEMS) may stick together and cease to work due to the attractive nature of van der Waals and Casimir forces. Recent 
^ years have also witnessed a huge activity in the development of metamaterials (MMs) [3], boosted by the possibility 
^ • that such engineered media ™^ give rise to novel optical properties at selected frequency ranges, including negative 
qh' refraction Q, perfect lensing and cloaking Q, among others. Generally speaking, metamaterials are made of 
I— I, micro and nanostructures carefully designed to collectively present a particular electromagnetic feature. They are 
able to produce the aforementioned striking phenomena, inaccessible with natural materials, due to the significant 
, magnetic activity built into them, starting at microwave frequencies and going all the way up to the optical range. 
■ Substantial magnetic activity at high frequencies is exactly what brings Casimir physics and metamaterials close 
' together. It is known for a long time that the attractive character of Casimir forces is not universal, that it actually 
00 depends crucially on the geometry and composition of a given setup. This naturally leads to speculations about 
the feasibility of customized systems displaying repusive Casimir forces, or maybe strongly reduced attraction, thus 
, providing a so desired anti- "stiction" effect. An interesting possibility of flipping attraction into repulsion was first 
' demonstrated by Boyer Q, where he showed that a perfect conducting plate repels a perfect magnetic one with 
On . vacuum in between. The requirement of having a perfect magnet was relaxed in subsequent works [9], where it 
' was shown that non-ideal materials could give rise to repulsion (provided one was mainly electric and the other 
^ mainly magnetic), but even this more realistic scenario is constrained by the absence of natural strong magnets at 
high frequencies t23i] . However, recent developments in nanofabrication have resulted in metamaterials with magnetic 
response in the visible range of the electromagnetic spectrum , fueling the hope for Casimir repulsion [ll[ . The 
^ . expectation is that, by tuning this magnetic response to the right frequency range and making it strong enough, one 
■ - ' could produce an experimentally measurable Casimir repulsion between, say, a MM slab and a thin metallic plate, or 
at least a significantly reduced attraction. 

The major issue about Casimir repulsion and, more generally, about any degree of Casimir force control is that 
the Casimir interaction between real dispersive materials is a broadband frequency phenomenon, as shown by the 
Lifshitz formula expressing the force between two semispaces as an integral over all (imaginary) frequencies with an 
exponential cut-off c/d [16]. In previous works [13, [3 '^^ have shown that this is a formidable obstacle on the route 
to Casimir repulsion, as typically the magnetic response of a MM is restricted to a narrow frequency band. In this 
work we would like to discuss an additional issue, briefly touched in (see also flQ*!), concerning the modeling of 
magnetic activity in a given metamaterial. As we shall see below, this is a very important discussion since the Casimir 
force depends crucially on the choice of the proper model. 
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II. METAMATERIALS: A SIMPLE MODEL 



As our first step, we want to present a simple yet effective procedure of how to characterize the magnetic activity 
of a metamaterial. This is better illustrated by tackling a specific example, so let us consider the metamaterial 
depicted in Fig. 1 It consists on a periodic array of double layered conducting cylindrical sheets, each one 

being characterized by its radius r, separation s between the sheets and length L. The layers themselves have some 
structure - they are tinily split apart on opposite sides, in such a way to prevent currents of completing an entire 
loop (see fig. 2). The array itself has period b, in addition we assume that the cylinders are long {r/L ^ 1) and that 
there is little space between the two layers (s/r <C 1) of a given cylinder. By turning on an external magnetic field 




^ y FIG. 2: Structure of each cylinder, seen 

from the top. 

FIG. 1: Sketch of the array of cylinders. 
The distance from center to center (the 
period) is b and each cylinder has radius r. 
For the sake of clarity, we don't show the 
structure of each cylinder in this figure. 

H{t) — H{t)z — iJoe~*'^*z parallel to the cylinders we induce circulating currents on them, that therefore induce a 
new magnetic field Hi [23] ■ The total field is therefore 



Ht(x) = H + H,(x) = {H + Hi{^))z , 



(1) 



where x is the position where the field is being evaluated and -ff/(x) is still to be determined. In this approximate 
model, our first assumption is that despite having different expressions inside and outside a given cylinder, is 
space independent in each of these two regions. Inside of a given cylinder, the field may be naturally divided into two 
contributions: the one that originates from that cylinder itself and another that comes from all the other cylinders 
combined, recalling here that the cylinders are finite in length and as such they "leak" some magnetic field to their 
surroundings. The first contribution is readily expressed in terms of the induced current per unit of length J with 
the help of Ampere's law, resulting in H^p = J/c. This of course holds exactly only for infinite cylinders, but is 
an excellent approximation for cylinders sufficiently long. For the second contribution wc shall present an heuristic 
argument, based on the conservation of the number of field lines. We know that the field strength in a given location 
is roughly proportional to the density of field lines in that location, so the first contribution gives 



H 



(1) 



kN 



(2) 



where N is the number of lines, pin is the density of lines and k is a proportionality constant. By following the same 
reasoning, in order to get the second contribution we must find the density of field lines "leaked" by all the cylinders, 
which is simply the density pout "leaked" by one cylinder times the number of cylinders A^cyi 



(3) 



where A is the (infinite) area of a plane perpendicular to the cylinders and we used the fact that a given (finite) 
cylinder leaks all of its lines to its surroundings. From ([3]) we get 



H 



Pi 

(1) 



7rr 



(4) 
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and, knowing that H'f^ and ijj^' have opposite signs, we may use our previous results to rewrite the total field inside 
a cylinder as 

1 1 Trr^ 

H'^^H+-J —J. (5) 

c c b'^ 

Since the field outside the cylinders may be worked out in a similar way, we shall just quote the result 

Hr = H--^J, (6) 

and point out that the only difference is that i?!^-* vanishes in this case. The key advantage of expressions (O and 
([S]) is that they are written in terms of the current density J, to what now we turn our attention. Each cylinder may 
be viewed as an RC circuit, and application of Faraday's law in one of them gives 



£ = 27rmJ - — - = — 

iLuL c at 



Ho H — J pj- J 

c c 0^ 



(7) 



where where a is the resistivity of the cylinder sheets and C is the effective circuit capacitance also per unit of length. 
By assuming that it evolves harmonically in time, i.e., that J = Joe"*"^*, we get 

'^"^ (8) 



27Tac c iuj_ / 1 7rr^\ ' 

Tzr iuCTTr^ c \ b'^ J 

which might be used in ([5]) and ^ to give explicit expressions to the fields. We are not, however, so interested in 
these absolute expressions as we are in their spatial average, since the latter presumably gives a better account of 
effective properties, specially at large scales. Starting from the integral form of Faraday's and Ampere's laws 

E dl = -i^ /b cIs , / H • dl = i / J • ds, (9) 



C Cdt Jg J (J C J g 

we see that the important quantities in these equations are the flux of B and the circulation of H. This clearly 
suggests that B is to be averaged over a surface while H should be averaged over a line. For a periodic system we 
need to average over an unit cell only, so, using the unit cell defined in Fig. 3 and the fact that Ht points in the 

Bavc - ^ ^ B • ds = 1 ^ Ht • ds = (10) 
Have = 7 ^ Ht • dl = H^^' , (11) 



b 



where Sz is either of the two faces perpendicular to the z-axis, and is one of the edges parallel to z. Now, following 
[20| and defining the effective magnetic permeability as 



we get finally 



Moff(c.') = (12) 

■Have 



where we defined of the filling factor /, the resonant frequency and the dissipation coefficient 7^. 

There are some peculiarities about the previous formula that we should point out. The first thing is the functional 
resemblance between expression (|13p and the Drude-Lorentz model for the electric permittivity, the only difference 
between them being the presence of a uP' in the numerator of (fT3|l . As we shall see in the next section, this leads to 
interesting consequences in the Casimir force. In addition, it is easy to see that /icff does not go to unity in the high 
frequency limit, as we should expect on physical grounds. This tells us that (|13p cannot be valid up to arbitrarily 
high frequencies, what makes perfect sense given that all our calculations are based on an effective medium picture 
that breaks down for very small wavelengths. Our final remark is that, despite the problems in the high frequency 
range, expression (jl3p is totally acceptable from causality requirements, since it is analytic in the upper half plane 
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FIG. 3: Illustration of the unit cell used 
in our array. Again, for the sake of sim- 
plicity, the cylinder is shown as a single 
sheet. 



of complex frequencies. Indeed, it may be shown that its real and imaginary part are linked together by a shghtly 
modified Kramers-Kronig formulae 

M'M = (l-/) + -Ppj/ /"i^^~\ 

Jo y - ^ 

M = P/ dy 14 

where — Re/i(aj) and fJ-"{Lu) = Im/i(w). 



III. THE CASIMIR FORCE 



Now that we have an expHcit formula for the magnetic activity of our MM, we may proceed with the evaluation of 
the Casimir force. There is, however, an important issue we have to deal with before going directly to calculations. 
The metamaterial described in the previous section is severely anisotropic, and ()13p gives the right magnetic response 
only when H is parallel to the z-axis. As the Casimir force depends on virtual fluctuations coming from all directions, 
the inevitable conclusion is that //off has to be generalized in order to take arbitrary propagation directions into 
account. Unfortunately, the simplicity present in the previous section disappears when we depart from the particular 
case where H//i, with effects like polarization mixing and spatial dispersion (2^ coming into play. In addition, this 
high degree of complication obscures all the physics we are trying to analyze, namely, magnetic effects affecting the 
Casimir force. So, in order to avoid unnecessary distractions, from now on we consider a toy model that is based on 
the one we have been discussing, but with an isotropic magnetic permeability described by eq. (|13p and an isotropic 
electric permittivity given by 

ecffM^l-^ , (15) 

where fl^, and 7e are respectively the electric oscillating strength, resonance frequency and dissipation coefficient. 

The Casimir pressure between two or more bodies can be calculated in fairly general situations using techniques 
such as the scattering formalism [l^], the Krein formula [l^, the argument principle [l3|, path integrals [IHl, and 
others. A thorough discussion of these formalisms is outside the scope of this paper, but for our purposes it suffices 
to say that the result obtained for the pressure is that it depends basically on the reflection properties of the bodies 
and the distances among them. For the particular case of two real (i.e., non- ideal) isotropic half-spaces separated by 
vacuum, as shown in Fig. 4, the pressure is given by 



where K3 — ^k"^^ + ^ jf? and the reflection coefficients i?™, i?™, are given by 



i?r(*^,k||) = \ , (17) 
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where e„ and /z„ are respectively the permittivity and permeability of the n-th half-space. Expression (|16p is slightly 
more general than the original Lifshitz formula |,16| , to which it reduces in the case of /i„ = 1 . 

In our conventions, a positive (negative) value of the force means attraction (repulsion). We are interested in the 
case where one of the half-spaces (see fig. 4) is purely metallic, implying in 

ei(u;) = l--^^^i , Mi(w) = l, (18) 

where Jli is the metallic plasma frequency and 71 is the dissipation coefficient, and the other half-space is a MM 
characterized by e2(w) = ecff('^) and H2{^) = Moff(w), given respectively by and 




FIG. 4: Illustration of the setup used in 
this paper to calculate the Casimir force. 



z = —d 



z = 



Before going straight to the results we would like to point out a small subtlety brought by the use of fi^s in the 
Lifshitz formula. Buried in expression pop is the assumption that the reflection coefficients are analytic in the whole 
upper plane of complex frequencies, which is based itself on the assumption that ei^2(*C): Mi,2(*C) ^ 0. It is easy to see 
that functions ei.2(*C) meet this condition for any choice of parameters, but for fJ-2{i£,) it only holds when r/b < I/a/tt- 
Fortunately, in our case we have necessarily r /b < 1/2 < l/\/7r so we don't have to worry, but this little analysis shows 
that one must be careful with the validity conditions of the Lifshitz formula when considering alternative permittivity 
and permeability models. 
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FIG. 5: The Casimir pressure between a 
metallic half-space and a MM modeled by 
permittivity ([15} and permeability ([T9)) . 
for several values of the magnetic oscillat- 
ing strength Qm ■ The parameters used are 
Qi/Q = 0.96, ji/Q. = 0.004, n^/Q = 0.04, 

0.005. The frequency scale Q. is chosen 
as 1.43 X lO^'' rad/s and A is defined as 
A = 2-Kc/Q.. 




d/A 

FIG. 6: The Casimir pressure between a 
metallic half-space and a MM modeled by 
permittivity (|15p and permeability (|13p . 
for several values of the filling factor /. 
The parameters used are the same as in 
figure 5. 



In Figs. 5 and 6 we compare the Casimir pressure as a function of the separation d for different MMs. We see that 
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when the magnetic part is given by a Drude-Lorentz contribution, namely 

/^dl(^) = 1 - , (19) 

the force decreases at least as fast as l/d'^, while for MMs modeled by we have a transition region where FjcS?' 
actually increases before it goes to zero for large distances. We see also that in the first scenario the force might 
even assume negative values for strong enough magnetic activity, while on the second case the enhancement of the 
magnetic part just makes attraction stronger and stronger in the transitional region. This is a clear indication that 
it is impossible to achieve repulsion with a MM described by (|13l) , and it is actually not difficult to check that this is 
indeed the case. As we can see from (jl6p . a necessary (but not sufficient) condition for a negative force to occur is 
that, for at least one polarization, i?j(z^) and i?2(*C) must have different signs. This information is entirely contained 
in the numerators of both reflection coefficients in ([TT]), and a simple (but tedious) analysis shows that when the 
permittivities and permeabilities of the materials involved are given by p^ . p8|) and lfT3|) we have 



^ni^m^ - + t^nm^nme/c^ > o , vfc||,^,n, (20) 

what prevents any difference in sign between R{ and i?2. This is to be compared with the case where the MM 
permeability is described by (fTO|) . since in this instance it is possible to have R'[^{i£^) ■ R2^{i£,) < and therefore 
Casimir repulsion is not ruled out. We see that a direct consequence of all this reasoning is that repulsion is guaranteed 
(at least for a given range of distances) provided /x(i^) is significantly larger than e{i^) and that attraction is enforced 
if <C e(iC)- This usually gives an useful rule of thumb to check if repulsion is possible, or if by changing the 

relevant permittivities and permeabilities we are approaching (or receding from) Casimir repulsion. 

That two similar models give widely different predictions for the force is a very interesting fact on itself, and, 
although we just discussed its mathematical roots, a physical explanation for it is still lacking. The problem is that 
the insertion of magnetic effects in the Casimir effect makes its physical analysis notoriously difficult, even the idealized 
Boyer setup Q is not easily grasped from a physical point of view. Part of the explanation for why materials described 
by either /Xeff(w) or /iDL(^^) are so differently affected by vacuum fluctuations is probably due to the fact that the 
former is a diamagnet for almost all frequencies (except around the resonance), while the latter is a paramagnet up 
to the vicinities of the resonance (and then becomes a diamagnet), but to go beyond that and state something more 
conclusive about their distinct physical behavior is, at the current stage, premature. 

We would like to end this section with a brief discussion about the effects of electromagnetic anisotropy. As it was 
discussed to some extent in [18], anisotropy usually works against repulsion, but in this particular case things may be 
a bit different. Let us consider an hypothetical MM that is electrically isotropic but has a tensorial magnetic response 
u described by (jl3p in one direction (say, the z direction) and is magnetically inert in the plane perpendicular to it 
[25j . This kind of MM gives rise to a slightly reduced Casimir attraction if compared to an isotropic MM described 
by (Uni), and the reason is that we have Ma;x(*C)/eoff(iO) Ai!/a(*?)/£cff(«C) > Mcff (iO/^off («6 (since ^ixx{iCj = fJ-yyi'iS.) = 
1 > /j,eff (j^C)); and hence, by the rule of thumb discussed earlier in this section, we should expect less attraction. 



IV. CONCLUSIONS 



In this paper we discussed how to obtain the magnetic permeability of a specific MM from a theoreti- 
cal/phenomenological point of view, and how the Casimir force is affected by it. Unfortunately, the kind of magnetic 
activity we obtained does not lead to Casimir repulsion, nor even to reduced attraction. Incidentally, these results 
also show that for Casimir purposes, extrapolations of models for electromagnetic properties must be done in a very 
careful fashion, as these extrapolations may be crucial to the determination of the Casimir force. Our final conclusion 
is that, as the description of practical MMs is typically known only on a narrowband of frequencies and incidence 
angles, much has yet to be learned before a definitive statement about the possibility of repulsive Casimir forces can 
be made. 
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